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ABSTRACT 

Context. Oxygen is the third most abundant element in the universe, but its chemistry in the interstellar medium is still not well 
understood. 

Aims. In order to critically examine the entire oxygen budget, we attempt here initially to estimate the abundance of atomic oxygen, 
O, in the only one region, where molecular oxygen, O2, has been detected to date. 

Methods. We analyse ISOC AM-C VF spectral image data toward p Oph A to derive the temperatures and column densities of H2 at the 
locations of ISO-LWS observations of two [OI] ^Vj lines. The intensity ratios of the (/=l-2) 63 /jm to (/=0-l) 145 /jm lines largely 
exceed ten, attesting to the fact that these lines are optically thin. This is confirmed by radiative transfer calculations, making these 
lines suitable for abundance determinations. For that purpose, we calculate line strengths and compare them to the LWS observations. 
Results. Excess [O I] emission is observed to be associated with the molecular outflow from VLA 1623. For this region, we determine 
the physical parameters, T and A^CHt), from the CAM observations and the gas density, n(H2), is determined from the flux ratio of the 
[Oi] 63 fim and [Oi] 145 /im lines. For the oxygen abundance, our analysis leads to essentially three possibilities: (1) Extended low 
density gas with standard ISM 0-abundance, (2) Compact high density gas with standard ISM 0-abundance and (3) Extended high 
density gas with reduced oxygen abundance, [0/H] ~ 2 x 10"^. 

Conclusions. As option (1) disregards valid [Oi] 145 /jm data, we do not find it very compelling; we favour option (3), as lower 
abundances are expected as a result of chemical cloud evolution, but we are not able to dismiss option (2) entirely. Observations at 
higher angular resolution than offered by the LWS are required to decide between these possibilities. 

Key words. ISM: abundances - ISM; molecules - ISM; dust, extinction - ISM; clouds - ISM; jets and outflows - ISM: individual 
objects: p Oph A, VLA 1623 



1. Introduction 

Oxygen is the most abundant of the astronomical metals and 
is, as such, of profound importance for the chemistry of the 
interstellar medium (ISM). Therefore, its role and its relative 
abundance in the various phases of the ISM should of course 
be known and understood. |Quan et al.| ( |2008 ', and references 
therein) present a recent overview for our understanding of 
the abundance of interstellar molecular oxygen. In general, the 
amount of O2 in molecular clouds has been below detection ca- 



pability of the dedicated space missions SWAS (e.g., Goldsmith 
et al.|2000t|Bergin et al.|2000 ) and Odin (e.g., Pagan i etal.|2003| 
Sandqvist et al. 2008 1. In merely one single location, viz. the 



dense molecular core A in the p Ophiuchi cloud, did Odin detect 
a weak O2 line at the ~ 5cr level ( |Larsson et al.|2007] l. 

The oxygen abundance in the interstellar medium has be en 
determined t o be about 3.0 x lO' "^ ( |Savage & Sembach|l996| to 
3.4 X 10 ( jOliveira et al.||2005| ). This is smaller than the solar 
value, viz. [O/HJg = 4.6 X lO""* ([Asplund, Grevesse & Sauval 



Send offprint requests to: R. Liseau 

* Based on observations with the CAM-CVF (Cesarsky et al. 
1996) and the LWS (Clegg et al. 1996) onboard the Infrared Space 
Observatory, ISO (Kessler et al. 1996). 



2005 [ [AUende Prieto"2008'p] indicating that in the interstellar 
medium (ISM), oxygen is depleted in the gas phase. 

Hollenbach et al.| ( |2009) l examine theoretically the oxygen 



chemistry in various interstellar cloud conditions in considerable 
detail, taking into account processes both in the gas phase and on 
grain surfaces. Of particular relevance in the context of this paper 
would be their results applying to weakly irradiated photon dom- 
inated regions (PDRs). Weak shocks also heat (and compress) 
the gas which, overcoming low-temperature reaction barriers, 
can potentially initiate a rather complex chemistry. Gusdorf et al. 
( 2008| l present results for oxygen bearing species, including O2, 
which are produced by grain sputtering in C-type shock models. 

The diagnostic usefulness of H2 lines for the study of 



shocked gas has been discussed in detail by [Neufeld et al. ( 2006 
and references therein). It had also earlier been shown that 
Herbig-Haro (HH) flows, which are optical manifestations of 
interstellar shock waves, are sttong emitters of [O i] 63 yum line 



emission (Liseau et al. 1997 1. In the present paper, we will com- 



Ayres et al. 1 2005 1 and 



Landi et al. 



1 2007 1 present evidence in 
favour of twice that value, i.e., 8.7 x lO""*. In a recent evaluation, 



Melendez & Asplund|( 200 8 1 concluded that a value of 5 x 10""* would be 
consistent with the existing different model atmospheres. There is gen- 
eral agreement, however, that the ratio of the carbon-to-oxygen abun- 
dance is unaffected by the different analysis techniques and remains 
close to one half. 
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bine observations in lines of these species of the star forming 
core p Oph A. 

There, one observes emission from both a PDR and from 
shocked gas and both processes will have to be considered a pri- 
ori. An overview of the region near the binary Class source 
VLA1623 ( ,Looney et al.||2000| ) in the pOphA cloud core, 
toward which spectral line data of both H2 and oxygen ex- 
ist, is shown in Fig.[T| This study focusses on the analysis of 
these data with the aim to obtain the oxygen abundance in the 
p Oph A cloud. 

This paper is organized as follows: In Sect. 2, we report 
on the ISOCAM observations, account for their reduction and 
briefly present the main results. In this section, we also give an 
account of the LWS observations and describe the results in con- 
siderable detail. In Sect. 3, we derive the temperature and the 
column density of the emitting H2 gas and present our results for 
the accompanying oxygen gas, with emphasis on the elemental 
abundance. Finally, in Sect. 4, we briefly summarize our main 
conclusions from this work. 




\ _ _ I 

16h26in45s 16h26m3(ys 16h26ni15s 16h2»m00s 

Right Ascension (J2000) 

Fig. 1. The orientation of the CAM-CVF frames in pOphA, 
with the pixel coordinate (0, 0) in the upper left corner, is 
shown superposed onto a partial Spitzer Space Telescope image 
(courtesy NASA/JPL-Caltech/Harvard-Smithsonian CfA, see 



Table 1. The peak position and the width of the Lorentzian pro- 
files used to fit the PAH features. 



PAH feature (yL^m) Peak position (cm ') FWHM(cm') 



6.2 
7.7 
8.6 
11.3 
12.7 



1596.9 
1295.8 
1176.5 
877.5 
781.2- 800 



70 
113 
40 
25 
40 



2. Observations, data reductions and results 

2.1. ISO-CAM CVF 

The pipeline processed ISOCAM-CVF data of pOphA were 
extracted from the ISO archive (TDT numbers 4520111 and 
45601809), covering a region of 3' x 3'. These observations trace 
the CO o utflow from VLA 1623 (e.g., ' Andre et al.|1990[[Dent et 
al. 1995| l in a map of 32x32 pixels, each representing 6"x6" (the 
angular resolution of ISO-CAM is diffraction limited). In each 
pixel of the image frames, the light has been dispersed by the 
Circular Variable Filter (CVF) to produce an infrared spectrum, 
so that the CAM data thus are represented by data cubes with 
two spatial and one spectral dimension. In addition, two wave- 
length channels result in frames, in which each pixel contains a 
short-wavelength (SW, 5 - 9.5yum) and a long-wavelength (LW, 
9.5 - 16.5//m) spectrum and the resolution is A/AA = 34-52 
( [Blommaert et al.|2003| . Each pixel had both of its spectra ex- 
amined individually and reduced manually as described below. 

Using known point sources in the region for positional ref- 
erence, it can be seen that the SW and LW frames are shifted 
by 3 pixels in the y-direction. With a 'dead row', in each of the 
S W and LW frames, there are ~ 850 spectra which cover the full 
5 - 16.5 yum range. The spectra of the two wavelength regions 
were joined together by adjusting either the S W or LW part ver- 
tically. The relative rotation was within a small fraction of a pixel 
and we did not attempt to correct for thaj^ 



2.1.1. Solid state and PAH features 

Each of these spectra shows presence of the silicate dust absorp- 
tion at lOyum, along with very strong emission at 6.2, 7.7, 11.3 
and 12.7 jim thought to be due to polycychc aromatic hydrocar- 
bons (PAHs, Leger & Puget|1984[[Ariamandola et al.|1985[ ). hi 
order to extract the emission due to PAHs, we fitted a continuum 



http : //gallery . spitzer . caltech.edu/Imagegallery/chrfefilpw the PAH features and over the silicate absorption (Fig.|2| 



php?cat=Astronomical_Images). The colours represent ef- 
fective wavelengths of 3.6 jum (blue), 4.5 fim (green) and 8.0 yum 
(red). The array size of 3' x 3' corresponds to 0.1 pc by 0.1 pc at 
the distance of the cloud (120pc: ,Lombardi et al.|[2008[ |Snow| 
et al.|200^ . A number of young objects are indicated, with the 
dominant B-star HD 147889 being off the figure, almost 15' 
westsouthwest of the center of the CAM-frame. The position of 
the Class source VLA 1623, which drives a collimated bipolar 
CO outflow ( | Andre et aL]|1990| ), is at the center of the dark 
flaring disk feature. The green colour, in particular, visualizes 
many flow phenomena, such as, e.g., HH 313 and the arcs of the 
southeastern flow. 



After subtracting the continuum, we fitted the PAH emission, as- 
suming that each line profile can be described by a Lorentzian 
profile ( Boulanger et al.|2005| l. The peak position and the width 
of each PAH feature are listed in Table[T] 

The star S 1 is illuminating the north-eastern part of p Oph A, 
while HD 147889 is heating the dust in the south-western region 
of the ISOCAM map (Fig.[T). The ratio of the PAH bands at 
7.7 yum and 12.7 fj.m can be used as an indicator of the degree of 



- More specifically, the equatorial J2000.0 coordinates for the SW- 
frame (TDT 45201 1 1) are RA = 16''26"21:'386, Dec = -24° 23' 59'.'28 
at spacecraft roll angle 98°31928. For the LW-frame (TDT 45601809) 
the corresponding data are RA = 16''26™21?259 and Dec = 
-24° 24' 11 "16 at roll = 98°5 1982. The coordinates refer to the respec- 
tive map centres. 
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Fig. 2. Left: The observed CVF-spectrum in CAM-pixel (3, 4), where PAH features dominate the spectrum. The short-wave and 
long-wave bands are shown as histograms with thin and thick lines, respectively. Fits to the continuum and the PAH features 
with Lorentzian profiles are shown in the upper two frames. In the lower frame, these components have been subtracted and the 
wavelengths of H2 lines are shown by vertical dashed lines for reference. Right: The spectrum of pixel (8, 22) corresponds to 
HH313 A and contains less dominant PAH-emission. Instead, a number of H2 lines are discernable. The upper two frames show 
again the fits to the continuum and PAHs, respectively, whereas the lower frame reveals the rotational line spectrum of molecular 
hydrogen. In addition to the emission features, there are also absorption bands from solids, viz. due to silicates at lOyum and due to 
CO2 at 15 //m. Theoretical models of both gaseous (red) and solid (black) CO2 absorption, at the appropriate spectral resolution of 
the CAM-CVF but arbitrarily shifted in the y-direction, are shown in the upper right corner of the middle frame. 



ionization of the PAH molecules and by comparing the maps of 
these two bands, it can be seen that the FUV photons from the 
early B-type star HD 147889 can penetrate the cloud and ionize 
the PAH molecules more effectively than those from S 1, a star 
of slightly later spectral type. This is entirely in line with the 
conclusions reached by Liseau et al. ( 1999 1 on the basis of their 
ISO-LWS observations. 

The data for the PAHs and the solid state features in p Oph A 
have been presented and discussed by Alexander et al. (2003 1. 
For the ice-ratio CO2/H2O, these authors give a range of - 0.4, 
but left unspecified which region in p Oph these numbers do re- 
fer to. In a limited number of pixels along the outflow, we de- 
tect a sharp absorption feature due to CO2 ice. At the resolu- 
tion of the ISOCAM-CVF, it is not possible to make a detailed 
study of the matrix of the ice, i.e., whether it is pure CO2 ice 
or mixed with water. However, we do not detect water-ice ab- 
sorption at 6.2 yum as any absorption would be filled in by the 



PAH 6.2 /urn emission feature. Also the broad libration band of 
water-ice would suffer from contamination by PAH emission at 
1 1.3 yum and 12.7 jjm. 



The broad absorption wing near 15/im could be indicative 
of CO2 gas and model fitting would require likely unrealistically 
large column densities (Fig.|2]i. This would be in large contrast 
to what has been observed elsewhere at much higher spectral 
resolution ( [van Dishoeck|[T998] ). Therefore, mainly due to the 
difficulty to correctly identify the continuum level we are unable 
to make any quantitative assessment of the column density of 
CO2 gas. The situation seems better for the ice feature though. 



for which we estimate A^(C02)i( 



(1.4 ± 0.4) X 10" cm- 



and where we have used the absorption cross section provided 
by Gerakines et al. ( 1995 1. This particular estimate refers to the 
pixel (8, 22), see Fig.|2j 
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Table 2. ISO-LWS [OI] line fluxes and line ratios for the VLA 1623 outflow and the neighbouring p Oph PDR 



TDT 


RA(J2000) 


Dec(J2000) 


10'^ F63^m 


S/N63 


10'^ Fi45^m 


S/N,45 


Project Name/ 


Comment 


number 


(h m s) 


(° ' ") 


(erg cm"^ s"') 




(erg cm~- s~') 




Position 




Flow 


















29200533" 


16 26 26.27 


-24 24 30 


18.52 ± 1.29 


14 


4.47 ± 1.50 


2.6 


VLA 1623 


Origin (flow 0" offset) 


29200534* 


16 26 17.22 


-24 23 06 


22.54 ± 0.80 


28 


4.22 ±0.86 


4.9 


VLA 1623 flow(l,l) 


150"NW(pa=304°) 


99900534- 


16 96 93 95 


-24 24 02 


18.52 + 1.27 


15 


4.47 + 1.35 


3.3 


VI A 1 693 flnwr9 1 ^ 


"iO" NW Tna— 304°1 


29200534 


16 26 29.29 


-24 24 58 


21.92 ± 1.00 


22 


8.94 ± 1.62 


5.5 


VLA 1623 flow(3,l) 


50" SE (pa=124°) 


PDR 


















45400801 


16 25 55.56 


-24 25 39 


10.08 ± 0.77 


13 


4.02 ± 0.80 


5.0 


ROPH_EW ew2 


~south 1' west 7' 


45400801 


16 26 08.74 


-24 25 40 


12.66 ± 0.62 


20 


4.92 ± 1.09 


4.5 


ROPH_EW ew3 


~south r west 4' 


45400801 


16 26 21.92 


-24 25 40 


10.56 ±0.81 


13 


3.45 ± 0.60 


5.8 


ROPH_EW ew4 


~south r west r 


45400801 


16 26 35.10 


-24 25 41 


12.66 ±0.36 


35 


4.41 ± 0.86 


5.1 


ROPH_EW ew5 


~south r east 2' 


PDRave' 






11.49 ±0.64 


18 


4.20 ±0.86 


4.9 






Offset' 


16 26 35.32 


-24 25 54 


9.88 ± 1.13 


9 


3.96 ± 0.93 


4.2 


VLA 1623 flow(4,l) 


150" SE (pa=124°) 


FlOW-PDRave 






11.05 




0.02 




552 = Line Ratio 


150"NW(pa=304°) 








7.03 




0.27 




26 


50" NW (pa=304°) 








7.96 




0.27 




29 


Origin (flow 0" offset) 


Flow-Offset 






12.66 




0.26 




50 = Line Ratio 


150"NW(pa=304°) 








8.64 




0.48 




18 


50" NW (pa=304°) 








8.64 




0.51 




17 


Origin (flow 0" offset) 


<FIOW-PDR,ve> 






8.68 




0.27 




> 25 




<FIow-Offset> 






9.98 




0.41 




24 





Notes to the Table: 

" Center position related to the strip map (TDT 29200534), 

' which is along position angle 124° (±180°) with 150" spacings of 4 individual pointings. 
^ The surrounding PDR provides one type of background estimate. 

'' TDT 29200534: This position of the strip scan is outside the flow, yet close enough to serve as another reference position for background 
estimation. 



2.1.2. H2-lines 

The residuals after fitting the PAHs show that in several pixels, 
there is a series of emission lines due to pure rotational transi- 
tions of molecular hydrogen, H2. The combination of SW and 
LW spectra made it possible for us to detect the transitions from 
S (2) up to S (7). These lines are not resolved at the resolution of 
the ISOCAM-CVF and were fitted using Gaussian profiles in or- 
der to estimate the line intensities (bottom right panel in Fig.|2|i. 
There are no detections of possible other line emitters, such as, 
for instance, [Ne II] 12.8 fim. 

The spatial distribution of the individual H2 lines is dis- 
played in Fig.|3] H2 line emission associated with the outflow 
from VLA 1623 is clearly seen. Along this flow, the emission is 
particularly prominent in the S (5) and S (6) lines. This indicates 
that the H2 outflow gas is in a state of elevated excitation and 
that temperatures are likely to be relatively high. In contrast, the 
PDR emission is especially prominent in the low-excitation S (2) 
emission, suggesting this gas to be at lower temperatures, consis- 
tent with the results from theoretical models of the p Oph A PDR 



([Liseau et al.|1999 ; Spaans & van Dis hoeck 2001 ;' Habart et al. 
12003 1 IKulesaeFal., 2005; Hollenbach et al.„2009; ). It should be 
noted that the sparse distribution of the S (4) fine is due mainly to 
the strong contamination by the 7.7 fim PAH feature. The bright 
pixels in the upper right correspond to the known Herbig-Haro 
objects HH 3 13 A and B (e.g., [Gomez et al.|2003 



2.2. ISO-LWS and [O I] lines 

The oxygen spectral line data analysed in this paper had been 
obtained with the Long Wavelength Spectrometer (LWS) on 
board ISO. In Table[2] the observations which were retrieved 



from the archive are identified by their TDT numbers. The pre- 
sented data are pipeline reduced and line fluxes were extracted 
from the detectors SW 3 for the [O i] 63 fim line and LW 3 for 
the [Oi] 145 jum line, respectively. To estimate the flux, we fitted 
the line after continuum subtraction with a Gaussian, the width 
of which was kept fixed to the resolution of the short- and long- 
wavelength LWS bands, i.e., 0.29 and 0.6 fim, respectively. The 
conservative error estimate reflects the peak-to-peak of the noise 
and the absolute accuracy of the LWS data is better than 30%. 
At these wavelengths, extended source corrections result in ef- 
fective LWS-beamsof 1.41 X 10"^ sr and 9.35 x 10"*^ ' 
lively ( [Giy et al.|2()03l ). 



sr, respec- 



The TDT 29200534 data set contains a strip scan along the 
position angle 124° and with 100" spacings. The center coor- 
dinates of the scan are given by TDT 29200533, i.e. toward 
VLA 1623. The pointings are along the CO outflow from that 
source. The data of the northwest flow provide the basis for 
our discussion below, whereas the observations toward the off- 
set 150" southeast provided a (quasi)simultaneous measurement 
of the background. This position, named "Offset" in Table|2]lies 
outside the CO outflow, which is changing its direction, yet it 
is sufficiently close to serve as a reference position for the flow 
data in the northwest. 

Although the use of a single reference position is common 
practice when observing dark clouds, we used also another set 
of observations, i.e. those of the surrounding PDR, to gauge the 
[OI] emission from and near the VLA 1623 outflow. The com- 
parison with this statistically estimated background should pro- 
vide us with the means to assess the accuracy of the results (cf. 
Table[3]l. These data showed that the PDR provides a background 
of essentially constant intensity and, consistent with the other 
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data set, that the 63 jum emission along the flow is enhanced by 
a factor of almost two, while the [Oi] 145 fim line resembles the 
background. 

For both data sets, the signal-to-noise ratio (S/N) of the 
63 yum data is generally much higher than 10, whereas that of the 
[Oi] 145 ;um lines is roughly 3 to 5. Subtracting the Offset and 
averaged PDR fluxes from the flow data revealed comparable 
[Oi] 63 fim excesses in both cases, but led to small residuals at 
145 fim, resulting in large relative line ratios, in excess of about 
20 (see the last two rows in Table[2]l. 

The [O i] 145 jum datum for the 50" southeast position shows 
an exceptionally large flux value. This cannot be explained in 
terms of contamination by the CO (7=18-17) line at 144.8 ;um, 
as no other CO emission is detected elsewhere in the spectrum 
(similarly, any potential blending of the [O i] 63 fim line with 
H20(8o8 - 7i7) can generally be dismissed on similar grounds). 
The continuum level is also enhanced, by a factor of 3 to 4, 
whereas the [O i] 63 yum emission is of the same order as the 
other flow values (Table|2]l. However, this position is outside the 
CAM frame containing the H2 map, and it is this map which is 
at the focus of the following discussion. 



3. Discussion 

3.1. Extinction of the H2 lines 

In dense clouds, corrections for extinction by dust might become 
necessary even in the infrared (Appendix A). The CVF admits 
the three para-Hnes S(2), S(4) and S(6) at 12.3, 8.0 and 6.1 fim, 
respectively, and the three ortho-lines S(3), S(5) and S(7) at the 
respective wavelengths of 9.7, 6.9 and 5.5 ;um. The integer in 
parentheses refers to the rotational quantum number of the lower 
state, 7io„ = Jup - 2. The location of the lines in the CAM-CVF 
spectra is shown in Fig.|2] 

The S(3) line is coinciding in wavelength with the broad 
lOyum-silicate feature and is therefore most sensitive to the dust 
extinction. In the observed CAM-CVF frame, both high and low 
values of Ay are derived (see Figs.|4]and|5]), with an average of 
21 ± 10 magnitudes over 63 pixels, in which we detected six H2 
lines. For a larger number of pixels, i.e., 162 pixels with 5 de- 
tected lines, this value is not significantly lower, i.e. 16± 1 1 mag. 
Partially responsible for the large spread might be the fact that 
the actual extinction curve is different from the one used here. 

Extinction values for extended regions of the pOph cloud 
and as high as those derived here have been reported also by oth- 
ers (e.g., [Davis et al .179991 [Whittet et al.|2008j l. Especia lly re- 
markable, and in line with our own findings, is the result of Davis 
let al.| (|1999') for HH313A, who found at 0'.'6 resolution varia- 
tions of the K-band extinction by 3 magnitudes (AAv=30mag) 
over only 6", i.e., the size of one CAM pixel. Obviously, the sur- 
face layers of the pOph A cloud are extremely inhomogeneous. 
In fact, this is evident already in the original of Fig. [T] which re- 
veals the large increase of the extinction even longward of 3 fim 
and also clearly shows the dust distribution over the flow region 
to be patchy and filamentary. 

The observed optical depth of the silicate absorption feature, 
Tsii, could potentially be used as an independent way of esti- 
mating the extinction (e.g., Whittet et al. 2008, and references 
therein). In the lower panel of Fig.|5] a plot of rsii versus Ay is 
shown. Evidently, no close correlation, such as that shown by 
the straight line and valid for diffuse clouds, does exist for the 
dense core p Oph A. From the work by Chiar et al.| (2007 ) it is 
evident that in dark clouds, this relationship does generally not 
hold for optical depths larger than about 0.6, corresponding to 



S(2) 




E ^. 




S(4) 



S(B) 






S(7) 



+ 



Fig. 3. For convenience, the ISOCAM-CVF frames have been 
rotated (cf Fig.[T]l and arranged so that the lines originating from 
the para-states are shown in the panels to the left and those from 
the ortho-states to the right. The dark horizontal bar is due to a 
row of reduced pixel sensitivity, whereas the dark edges are due 
to loss of data by the alignment procedure. The image dimen- 
sions are 192" x 210" in the horizontal and vertical direction, 
respectively. In the lower right corner, the PDR is clearly seen, 
especially prominently in the S (2) line. The scale bar is in units 
of erg cm"^ s ' and the brightest spot is HH313 A, the proper 
motion of which associates it with the CO outflow. To its right is 
HH 313 B, which is likely unrelated to VLA 1623 but rather em- 
anates from GSS 30 ( Caratti o Garatti et al.|2006j ). The position 
of the outflow source VLA 1623 is indicated by a white cross. 



Ay ~ 10 mag. In our work, Tsii generally exceeds 0.6 (Fig.[5]l and 
we decided against the use of the diffuse cloud relation. 

3.2. Pliysical parameters of tlie H2 gas 

The H2 parameters presented here are based on a 'rotation dia- 
gram' analysis, in which it is assumed that the emitting gas is in 
local thermal equilibrium (LTE) at a single gas kinetic temper- 
ature and that the emission is optically thin (line center opacity 
much less than unity). In Appendix B, the method and the valid- 
ity of the assumptions are examined. 



6 



R. Liseau and K. Justtanont: Oxygen in dense interstellar gas 



\ 

CM 



10' 



101 



101' 



101' 



- 

1 

' Hg-line S(2) 


1 

S(3) 


1 

S(4) 


S(5) 


1 

S{6) 


1 1 

S(7) 


- 
- 














- 

- 














: 


- 












- 


r 




<> 
















* 














































m 




_^(^6m) 12.3 


9.7 


8.0 


6,9 


6.1 


5.5 




: (mag) 1.1 


2.3 


0.9 

. 1 


0.9 


1.0 

. 1 


1.0 





2000 4000 6000 8000 

Upper Level Energy, E, /k (K) 

up 

Fig. 4. The rotational H2-diagram for CAM pixel (8, 23), 
HH 3 13 A. At their respective upper level energy (in K), the tran- 
sitions are identified. Points with error bars (para: red diamonds, 
ortho: blue squares) refer to observed (open symbols) and ex- 
tinction coiTected values (filled), respectively. The value of the 
visual extinction, derived from -minimization, is Ay = 28 mag 
and the corresponding A i in magnitudes are given below the 
graph. The inverse slope of the straight line is a measure of the 
gas temperature, T = 1050 ± 150 K. At this particular location 
in the outflow from VLA 1623, the derived ortho-to-para ratio is 
o/p = 2.1+0.2. 



The H2 gas associated with the VLA 1623 outflow is remark- 
ably isothermal, T = 10^ K (formally 962 + 129 K, cf. top panel 
of Fig.|5]l and the average column density of the warm gas is 
A^(H2) = 3.5 X 10'^ cm"^. The ortho-to-para ratio along the flow 
is determined as 2.2 ± 0.4, respectively. 

The mass of warm H2-gas along the outflow from VLA 1623 
amounts to > 8 x lO"'* Mq, which radiates more than 7 x 10"^ Lq. 
This mass is nearly comparable (~ 50%) to that of the cold gas 



traced in CO in the northwestern outflow as derived by Andre et 



al. ( 1990|), assuming a 25% smaller distance ( 1 20 pc : ,Lombardi 
et al.|2008)|Snow et al.|2008] ). These estimates exclude the con- 
fusing blueshifted flow from SM IN ( Narayanan & Logan 2006 1. 



3.3. CP excitation and predicted line fluxes 

The excess flux in the [O i] 63 fim line is limited in extent to 
the VLA 1623 outflow, delineated by the rather uniformly dis- 
tributed rotational H2 line emission. Especially at the LWS- 
position 150" NW are both H2 and [OI] emission directly seen 
to coincide spatially (Fig.[6]l. Inside the cloud, the atomic oxygen 
is neutral and the warm gas seen in the rotational H2 lines is ex- 
pected therefore to contribute significantly to the [O I] emission 
detected by theLWS. 

Gas at temperatures around 10^ K can be expected to emit 
prolifically in oxygen fine structure lines (see Fig.|7]i. Toward 
150" NW, the observed flux ratio 7^63 /^145 » 10, which in- 
dicates that the emission is optically thin. However, we did in 
general not assume optically thin emission, but treated the ra- 




20 40 
Visual Extinction, (mag) 

Fig. 5. Top: Visual extinction. Ay, and temperature of the H2- 
gas, showing no dependence as expected. Bottom: Visual ex- 
tinction. Ay, and optical depth in the silicate absorption feature, 
Tsii. The error symbol shows the average of the points and their 
standard deviation. The dashed line displays the relation for dif- 
fuse clouds of |Whittet et aL] ( |2008l ). 



T(K) 




Fig. 6. The spatial distribution of the kinetic gas temperature de- 
rived from the H2 S (2) to S (6) lines. The position of the outflow 
source VLA 1623 is indicated by a white cross and the LWS 
pointings are shown by circular beams of diameter 70" and their 
positions refer to 0", 50" NW and 150" NW (see Table|2}. 



diative transfer properly in the LVG approximation (for details, 
see' Liseau et al.|2006| l. We adopted the H2-collisional rates from 
Jaquet et al. } 1992 1 also for these computations and determined 
consistently the rates for the average ortho-to-para ratio of H2 
and accounting for 20% of neutral helium. 
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Table 3. Predicted [OI] fluxes for the ISO-LWS at position 150" NW, assuming [O/H] = 3.4 x 10 1 
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T 
(K) 


(cm-2) 


mo) 

(cm-2) 


M(H2) 

(cm-') 


^63 /.im 


(erg cm"2 s"') 




f 145//m 

(erg cm-2 s-') 


F 63 IF 

pred" 


^obs / ^pred 

63yum'^ 


^obs / -^pred 

145 pm* 


20 


3.0 X 10^2 


2.0 X 10" 


3x 10' 


110 


4.2 X 10-'* 


9.7 X 10-" 


4.5 X 10-'« 


92 


262 


44 


1000 


3.5 X 10" 


2.4 X 10'* 


5 X 10' 


0.12 


8.9x 10-'2 


-0.06 


8.0 X 10-" 


11 


1.1 


0.025 








3x 10' 


0.05 


1.5 X 10-'" 


-0.01 


4.5 X 10-'2 


33 


0.07 


0.004 



Notes to the Table: " Observed ratio Fas/^ws ~ 50 (see Table|2]l. * Observed fluxes Fes = 1 x IQ-" and F,45 = 2 x lO-'-* erg cm-^ s-'. 




Temperature, T (K) 



Fig. 7. The predicted flux into the LWS beam of [O i] 63 /im 
(sohd Hnes) and [O i] 145 jiva (long dashes) is shown as function 
of the temperature and for the indicated values of the H2 density, 
assuming an oxygen abundance of 3.4x 10-"* relative to H and the 
column density determined from the rotational H2 observations, 
i.e. A^(H2) = 3.5 x 10'''cm-2. Values observed at the VLA 1623 
outflow position 150" NW are found inside the dashed lines and 
bound by the ±2cr values of the [O i] 63 jim line flux and of the 
temperature. The upper level energies for the [O i] 63 pim and the 
[O i] 145 //m lines are 228 K and 326 K, respectively. 



For the analysis of the oxygen lines, we used the values 
of T and iV(H2) determined from the observation of the H2 
lines. Regarding the volume density, n(H2), values in excess of 
10^ cm-^ are required by the very large flux ratio of the 63 /im 
to 145 //m lines, i.e. n(H2) ~ 3 x 10^ cm-"* (Fig.js]). This value 
refers actually to the warm gas, but it is also identical to that de- 



termined by Johnstone et al. ( 2000 1 from submm-observations 



of the cold dust. However, in the oxygen line analysis below, 
we examined a range in densities (and temperatures), whereas 
the oxygen abundance w as kept constant at i ts ISM value, i.e. 
3.4 X 10-"^ per H-nucleus ([Oliveira et al.|2005|). 



Fig. 8. The line flux ratios of the [O i] 63 //m and [O i] 145 //m 
lines are shown as a function of the temperature for five values 
of the H2 density. The slightly different beam sizes of the LWS 
have been taken into account and the computations assumed 
the common line width of 1 kms-'. Values observed along the 
VLA 1623 outflow are found inside the dashed lines and bound 
by the +2cr values (+300 K) of the temperature. 



3.3.1 . Cold gas phase oxygen 

A priori, one might not wish to dismiss the possibility that some 
or all of the [OI] line emission originates in the dense and cold 
gas, possibly associated with the high velocity CO gas, and we 
examine this case in the present section. 

As a conservative upper limit to the column density of the 
cold gas we take that for the corresponding average Ay, i.e. 
A^(H2) < 3 X 10^^ cm-^. This value is consistent with the find- 
ings also by others ( [Larsson et al.|2007l and references therein). 
Assuming that the cloud depth is comparable to the width of the 
CO outflow, i.e. 50" or 9 x 10'^ cm, the volume density is, again, 
«(H2) = 3 X 10^ cm-\ 

The upper level energies for the [O i] 63 /^m and the 
[O i] 145 /^m lines are 228 K and 326 K above the ground. As also 
shown in Table[3] oxygen gas at low temperatures, say T = 20 K, 
and with column density A^(0) < 2 x 10'^ cm-^ can therefore 
not account for the observed emission, in particular not for that 
of the higher excitation [Oi] 145 //m line. As can be seen from 
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this table, the 63 fim Une is very optically thick, 



10^ 



which would result in highly uncertain estimates of the oxy- 
gen abundance. In contrast, the 145 jim line is optically thin, 
Tus^jtn - 10""', and the theoretical flux in that line can be ex- 
pected to be more reliable. However, even for this large column 
density, the predicted line strength is inconsistent with the ob- 
servation by large factors. Lowering the average density would 
increase these inconsistencies. In addition, for the densest part of 
the pOph Acore, Liseau et al. (2006 ) estimated an atomic oxy- 
gen column density of only lO'' cm"^, decreasing the predicted 
line flux by another two orders of magnitude. In summary, cold 
gas emission, associated with that seen in low-lying rotational 
transitions of CO, does not contribute to the observed [O I] fluxes 
at any level of significance. 



3.3.2. Warm gas phase oxygen 

For the warm gas, we used the H2 parameters T = 10^ K and 
N(ii2) - 3.5 X 10''^ cm"-, assuming an oxygen abundance of 
3.4 X 10"^ relative to H. As a free parameter, the H2 density was 
varied in the range lO-' to 10* cm"-'. We assumed the common 
linewidth of 1 km s"' for both lines (FWHM). This corresponds 
essentially to the thermal width at lO-' K, providing maximum 
opacity in the lines. If associated with the outflow, their widths 
are likely broader 

In Figs.|7]and|8] the dependencies of the [O I] line fluxes and 
their ratios on the temperature and density are shown. Displayed 
are the rather wide ranges of 300 - 3000 K and 5 x 10"* to 
10* cm"^, respectively. Evidently, above the upper level energy 
of the [Oi] 145 yum line of 326 K, the dependence of the line 
fluxes on the temperature becomes decreasingly weaker. For ex- 
ample, allowing for variations with, e.g., the standard devia- 
tion of 150 K about the average temperature of 10^ K leads to 
changes in both lines by no more than 10%. Therefore, temper- 
ature variations are unlikely to significantly alter the results of 
Tabled 

The line center optical depths are also presented in this table, 
demonstrating that the lines are indeed optically thin. Therefore, 
any signifcantly different value of the column density would di- 
rectly alter the fluxes correspondingly. Such changes of A^(H2) 
would, however, be difficult to be reconciled with the H2 data. 

Also indicated in Fig.|7]is the parameter space occupied by 
the observed [O i] 63 yum line at the LWS position 150" NW and 
which corresponds to the ±2cr values of both the flux and the 
temperature. As it turns out, a unique solution is difficult to ob- 
tain and the results for two particular solutions are also presented 
in Table[2] In the following, these will be examined in more de- 
tail. 



3.4. Extended low-density gas and standard 0-abundance 

The first of these concerns gas at the rather low densities of be- 
low lO'^cm"^ and which is identified inside the dashed box of 
Fig.|7] This [Oi] 63 //m emission, filling the LWS beam, would 
also be in nice agreement with the observed H2 distribution. 
Sampled with 6" pixels, this H2 emission is clearly spatially 
resolved, showing some large-scale and coherent structure at 
~ lO-* K. In this scenario, the derived oxygen abundance would 
be essentially the one, which has been assumed by the calcula- 
tions, viz. [O/H] - 3.4x 10""*. Somewhat puzzling, though, is the 
non-detection of the accompanying 145 yum line at the level of 
10"'^ erg cm"^ s"' (see Fig. 7|. As a consequence, the observed 



flux ratio, F^/F 145, is much larger than what can be accommo- 
dated by any low density gas, radiating at the observed flux level. 

Given the evidence provided by the other outflow positions 
observed with the LWS, it seems not justified to simply dismiss 
the validity of the 145 fim data (see Table|2|i. The S/N of these 
lines is not overwhelming, which could cast doubt on the good- 
ness of the background corrected data. In general, when sub- 
tracting two nearly equal numbers, this creates a problem for 
the reliability of the result. However, in our particular case, sys- 
tematically higher [Oi] 145 yum flux levels (5.5 x 10"'^ erg cm"^ 
s"', rather than 4.5 x 10"'^ erg cm"^ s"', prior to the background 
subtraction) should become easily detectable. 

3.5. Compact high-density gas and standard 0-abundance 

A second possibility concerns gas which actually fulfills the line- 
ratio requirement, i.e. gas which would have to be at densities 
exceeding 10^ cm""* (Fig.[8]l. The observation of similar ratios 
also at the other positions along the flow attests to the likely re- 
ality of this result. Comparing to the data of Table|3] a beam 
filling below the 10% percent level would be indicated, which 
would correspond to linear source sizes of the order of 20". This 
is not unreasonable, as small scale structure is observed in ro- 
vibrationally excited H2 emissi on ([Pent et al.||1995] Davis et| 
al. 1999; |Caratti o Garatti et al.||2006|l, and also in [SH] Unes 
from HH 313 ( |Wilking et al.|1997||G6mez et al.|2003| [Phelps & 
Barsony 2 004| ). 

These structures are small inhomogeneities (^ 10") in the 
tenuous outflow gas, which are already pre-existent or have been 
broken off the cavity walls by Rayleigh-Taylor type instabili- 
ties. These experience impacts of essentially normal incidence 
and are shock-heated to several 10^* K. The non-detection of 
the [Nell] 12.8 /im line at 6" resolution (Sect. 2 .1.2) limits the 
shock velocities to below 60kms"', however (Hollenbach & 



McKee 1989). These HH shocks have been modeled in detail 
by Eisloffel et al.'(2000). Based on the observations in the NIR 
and in the optical, we know that their filling factors of the LWS 
beam are extremely small, i.e. 10""* - 10"^. We cannot exclude 
the possibility that these regions are larger, but hidden by extinc- 
tion. 

For these compact sources of emission, derived temperatures 
range up to some thousand degrees (1"- 10" scale), higher than 
the 10^ K obtained in this paper (50"- 100" scale). However, 
theoretical fluxes would not change by much, as the temperature 
dependence remains essentially flat. For instance, the increase in 
[O I] line emission from high density gas at about 2000 K would 
be below the 20% level. The effects from a reasonable range of 
temperatures are well within the 2cr flux boundaries considered 
in Fig. I?! and would, as such, not significantly alter the results of 
TablellT 

3.6. Extended high-density gas and low 0-abundance 

A caveat of this latter scenario is the observed extent of the H2 
emission as shown in Fig.|6] Line emission from accompanying 
atomic oxygen gas would fill at least half of the LWS beam, i.e. 
more than what is indicated by the /^obs/^pied ^ 0.1 in Table[3] 
This could potentially indicate a reduced oxygen abundance, i.e. 
formally to 2.4 x 10"^, corresponding to a column density of 
A^(0) = 2 X 10'^ cm"2. 

In this context, one issue concerns also the very origin of 
the atomic oxygen in this region. This gas is clearly associ- 
ated with the CO outflow from VLA 1623. This flow appears. 
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at least partially, to be embedded in the p Oph A cloud, the sur- 
face layers of which have been modeled as a PDR. According 
to the recent model by HoUenbach et al. ( 2009]) of the pOph 
PDR, the O" abundance is strongly reduced from its initial value 
at akeady quite modest values of Ay into the cloud. These au- 
thors model successfully the Odin results for molecular oxygen, 
O2, in pOphA ( Larsson et al.|[2007) l, the O2 being at its peak 
abundance. For their specific model and standard value^in their 
Eq. (32), we estimate an X{0) - 2.4x 10"^. This particular value 
refers to a visual extinction Ay ~ 4 mag and for higher values, the 
oxygen abundance is predicted to decrease further In addition, 
its coincidence with the value of the previous paragraph should 
be viewed as accidental, as the densities differ by a factor of six. 

To summarize, it is conceivable that the atomic oxygen ob- 
served in the outflow could be in dense cloud material with 
a reduced abundance compared to its initial ISM value. The 
emission would result from highly oblique shocks, as the flow 
is sliding along the cavity walls with shock velocities below 
5kms~\ resulting in temperatures of the dense wall gas up to 
10^ K. Hardly any results from theoretical models for shock ve- 
locities below 5kms"' are available (e.g.,iKaufman & Neufeld] 
1996] l. The extended emission in pure-rotational H2 lines origi- 



nates from this gas. This warm high density gas also dominates 
the [0 1] emission. 

3. 7. Summary: the oxygen abundance in p Oph A 

In conclusion, on the basis of the available observational ev- 
idence it is not possible to arrive at a unique solution to the 
oxygen abundance problem. If one disregards the [Oi] 145 yum 
data as too inaccurate, and therefore as irrelevant, the combined 
H2 and [OI] observations lead to an abundance which would 
be consistent with standard (depleted) ISM values. If, on the 
other hand, the [Oi] 145 //m data are invoked as relevant, the 
combined observations can be interpreted as indicating a con- 
siderably lower X{0). The latter scenario would be qualitatively 
in accord with recent models of cloud-FDR chemical evolution. 
To decide between these possibilities one would need to resolve 
the [OI] emission regions on scales smaller than 50" and, ide- 
ally, at the resolution of the H2 data with the ISO-CAM, i.e. 6". 
At 63 fim, the R4CS instrument aboard the upcoming Herschel 
Space Observatory with its 9" pixels should come close enough 
and should be capable of providing this information. 

4. Conclusions 

Below, we briefly summarize our main conclusions, which are 
based on the interpretation of observations obtained with instru- 
ments aboard ISO. 

• The analysis of ISOCAM-CVF specti-a of the VLA 1623 out- 
flow revealed that H2 is appreciably excited mainly along the 
flow. 

• The emission in the six rotational Unes S (2) to S (7) arises in 
gas of relatively constant temperature and column density. 

• Observations with the LWS along the flow and neighbour- 
ing regions show enhanced emission in the fine structure line 
[O i] 63 fim in the warm outflow gas. 



In pOph, the ortho-H20 abundance (relative to the number of 
hydrogen nuclei) in the gas phase h as variously been determined as 



< 3.8x 10 ^ 2.2 X 10"^ and 6.5 x 10"' (Liseau & Olofsson 1999 



Ashbyl 



|et al. |2000{|Franklin et al.|2008| respectively). These values refer to an 
average over a large area and/or a large beam and it is conceivable that, 
along the outflow, H2O abundances could locally be higher. 



Invoking a model for the [O I] line excitation and radiative 
transfer leads to estimates of the atomic oxygen abundance 
in the dense core p Oph A. These calculations use primarily 
the parameters determined at 6" resolution for the H2 gas, 
but examine also a considerably wider range in temperature 
and density. 

This results in three main scenarios: (1) low density extended 
gas (> 70") with standard ISM oxygen abundance; (2) high 
density compact sources (< 20") with standard ISM oxygen 
abundance; (3) high density gas (> 50") with largely reduced 
oxygen abundance. 

The third option would be qualitatively in accord with recent 
models of chemical PDR-cloud evolution. 
Observations with the Herschel Space Observatory can be 
expected to provide a test in the near future. 
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Appendix A: The extinction curve 

A widely exploited extinction curve is the 'average galactic law', calibrated for 
the near infrared by, e.g., Rieke & Lebofsky 1 1985 1. However, longward of about 
5 fim, no such generally accepted curve exists (e.g., Weinga rtner & Drainej200l"} . 
Dust properties may vaiy widely from location to location and unique extinction 
curves may be difficult to obtain. 

The optical depth of the extinguishing dust, at the wavelength A and along 
the line of sight, is given by r,i = 0.4A,i/loge, where the extinction at any 
wavelength, Aj, is given in magnitudes. The con'esponding extinction in the V- 
band. Ay, can be found from = 4.21 X 10^^ k,{A\j, where is the mass 
extinction coefficient (absorption + scattering) in cm^ g"' . For we adopt the 
values of |Ossenkopf & Henning| ( 19 941 bare grains, MRN, n=0) and normalize 
the extinction in the K-band (2.2//m) through Ak = 0.122Av, assuming a gas- 
to-dust mass ratio of one hundred and the column density-Ay relation of |Bohlin| 
etal. (1978l,i.e.,A'(HI) + 2A'(H2) = l.SxlO^'Av (in cm^^ mag"'). Henceforth, 
we will refer to this extinction curve as the OHBG curve. 

This coefficient for Ak for the OHBG curve is close to that given by 
|Rieke & Lebofsl^ jl985| i.e., 0.1 12, see Fig.|A.l|. It agrees excellently with 
the value suggested by [Savage & Mathis 1 1 979 viz. 0.123) specifically for 
the pOph cloud. The IRieke & Lebofsky (19851-curve is consistent with the 
A'(H)-A,i calibration in the I-band (0.9 ftm) suggested by Cardelli et al. 1 1989 1 to 
be applicable also for dense clouds (see also, ^Kenyon et al.|1998^ and provides, 
as such, a reasonable extrapolation into the infrared. 

The extinction in the pOph cloud is known to be 'non-average', with a 
larger ratio of selective-to-total extinction i?v ='^\/Eb-w ~ 5. The anomaly of 
the p Oph extinction in the optical and the UV has commonly been explained as 
being the result of a preponderance of larger grains. This could be due to grain 
growth produced by coating of small grains with icy mantles. The low gas-phase 
H2O abundance observed by SWAS (Ashby et al. 20001 would be in support 
for such a scenario and one might expect to see the growth effect by comparing 
observations with theoretical extinction curves for ice coated grains. However, 
it may be surprising to find that the observed width of the silicate absorption in 
pOphA, i.e., FWHM(10/im)=1.4/jm, is essentially that of 'average' dust, i.e., 
significant broadening due to ice opacity is not observed (see Fig. |A. H . 

Whereas dift'erences between individual extinction curves can become huge 
in the ultraviolet, errors introduced by a particular choice of curve are expected 
to become smaller in the infrared (Cardelli et al. 19891. 

For the OHBG curve, the extinction in units of magnitudes of Ay was varied 
by steps of 0. 1 mag. All H2 lines in a given location (CAM-pixel) were assumed 
to have suffered the same amount of extinction. Linear regression fits to the H2 
line flux data, with their error estimates, was made in a rotation diagram, re- 
sulting in a quantitative estimate of the goodness-of-fit, expressed equivalently 
by a^^ -value (Press et al. 19861. The best fit was then selected as that with the 
smallest providing the corresponding Ay for that pixel (see Fig.|4j. 



Appendix B: I-I2 excitation 

B.1. The single-temperature approximation 

The Herbig-Haro objects HH 3 1 3 A and B (Fig.[TJ signal the presence of shocks 
with strong temperature gradients on relatively small spatial scales jPavis et al.| 
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Fig.A.l. The extinction curve of Ossenkopf & Henning 



representing bare grains (OHBG) is shown as filled blac 



(1994: 



c dots 



and the interpolating solid line. The re-scaled data of Rieke & 
[Lebofsky (1985 ), using /c^ = (Aa/Av)/3.88 x 10"^ are shown as 
blue squares. Also shown are the models for ice-coated grains 
(green dots: thin ice mantles; red dashes: thick ice mantles). The 
wavelengths of the H2 lines admitted by the CAM-CVF are indi- 
cated by vertical bars. Toward the outflow from VLA 1623, the 
observed absorption feature near IQfiva has an average width of 
1.4jum (FWHM), i.e., essentiaUy that of the OHBG extinction 
curve. As is evident from the other curves, ice coatings would 
broaden that feature significantly. 



[T999l|Eisl6ffel et al.|200^ . The extent of these locahzed regions is comparable 
to (or smaller than) the size of the 6" CAM-CVF pixel, i.e. lO'^cm. In a ro- 
tation diagram, the data would not follow straight lines but display significant 
curvature, i.e., upturns implied by the higher temperature gas. 

The fit of the H2 line data to a straight line would ignore any such curva- 
ture, which could be expected for a distribution in temperature within the pixel 
field of view. Significant amounts of gas at higher temperatures would result in 
significant flattening toward the higher upper-level energies and would result in 
lower quality fits (cf. Fig|4j. This is not observed, however, when truncating the 
S (7) line from the data set with detections of all 6 H2 lines. The comparison of 
the 6-line and 5-line data sets in the region of overlap (63 pixels) demonstrates 
that the temperatures are identical within the observational errors. The overall 
mean temperature ratio is 1.00 ± 0.07. 

Selecting the data set with only the S (2) to S (6) detections nearly triples 
the number of pixels from 63 to 162, significantly enlarging the region of H2 
emission which can be analysed, i.e. Vl62 X 6", which is comparable to the 
angle subtended by the LWS (~70", [Gry et a l. 20031. Remarkably, also over 
this extended scale, the gas temperatures exhibit a uniform distribution: whereas, 
e.g., the extinction varies widely between a few and several tens of magnitudes of 
Ay, the temperatures stay rather flat at ~ 10^^ K (Fig.js}. Hence, we conclude that 
the extended warm H2 outflow gas near VLA 1623 is reasonably well represented 
by the single temperature approximation (see also Fig.|6}. 



|et al.|1999} . Data are available for collisions with HI, He, ortho-H2 and para- 
H2. The total collisional rate at a given temperature is therefore the sum of the 
fractional contributions by the collision partners. 

pOphA is a dense core, so that densities can be expected to be high and 
that the assumption of LTE therefore would be plausible, at least for the lower 
/-transitions (Fig. |B.l^ . Generally, in regions, where temperatures would be rel- 
atively low (T ~ 500 K), the critical density of the S (7) line would exceed 
10* cm"' and the line would likely not be thermalized. In contrast, at the nominal 
temperature of lO' K in pOphA, ncrit for the S (7) line would be below the av- 
erage density of the VLA 1623 region, i.e., n(H2) = 3 X lO' cm"' (see Sects. 3.5 
& 3.6) and we conclude that the assumption of LTE is justified. 
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Fig. B.l. The critical density, Wcrit = ^ui/Tui(r), for rotational 
transitions in the vibrational ground state, {v" , v' - 0, 0), for 
three values of the temperature, T . Considered are collisions 
with H I, He, ortho-Ha and para-Ha. The two extreme cases, i.e., 
for o/p - 3 (red and blue symbols) and o/p = (green dashes), re- 
spectively, are shown. Connecting lines have no physical mean- 
ing and are inserted for clarity only. In addition to the ortho-to- 
para ratio, parameters are X(H2) = 0.5 (all hydrogen molecular) 
and X(He) = 0.1. 



B.3. The optical depth of the H2-lines 

It is easy to show that the spectral lines are optically thin, t{v()) <s: 1, as long as 
the column density of warm H2 gas obeys ^(Ht) <*: 1.5 x 10~^ cm"^. 



The numerical coefficient 



based on the following 



max(n„/«tot) = 1-0, ininAii = l.Okms"', minT = hvo/k » 500K for the S(0) 
line and maxX(H2) = max A'(H2)/W(H) = 0.5. Realistic values for higher tran- 
sitions will always result in lower column densities, hence smaller optical depths 
(7" > 500 K) and, for the proper parameters, this condition is always checked a 
posteriori. Also, in outflow gas, Av » 1 kms"', decreasing t{vo) per unit col- 
umn even further The condition of optically thin emission is clearly fulfilled. 



B.2. Thermal equilibrium of the Hz-molecules 

The level populations become thermalized at densities which are higher than 
the 'critical density', Jicrit = Aii\/y\i[(T), where radiative and collisional de- 
excitation rates ar e equal (Fig.|B.l). In tha t figure, Einstein-A values have 
been adopted from [Wolniewicz et al.|{r998) and the rate coefficients for col- 
lisional de-excitation, yul(^)j were taken from the web-site maintained by 
D. Flower, http: //ccp7 . dur . ac .uk/ccp7/cooling_by_h2/ (see |Le Bourlot| 
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